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Continuous variable entanglement between two modes of a radiation field is usually studied at 
optical frequencies. As an important step towards the observation of entanglement between propa- 
gating microwave photons we demonstrate the experimental state reconstruction of two field modes 
in the microwave domain. In particular, we generate two-mode correlated states with a Josephson 
parametric amplifier and detect all four quadrature components simultaneously in a two-channel 
heterodyne setup using amplitude detectors. Analyzing two-dimensional phase space histograms for 
all possible pairs of quadratures allows us to determine the full covariance matrix and reconstruct 
the four-dimensional Wigner function. We demonstrate strong correlations between the quadrature 
amplitude noise in the two modes. Under ideal conditions two-mode squeezing below the standard 
quantum limit should be observable in future experiments. 



State tomography for more than a single mode of a 
radiation field allows to characterize photon sources that 
display entanglement between propagating photons. At 
optical frequencies well-established multi-mode state re- 
construction techniques based on single photon counters 
exist. These have already allowed for a variety of exper- 
iments that demonstrated entanglement and the EPR 
paradox for continuous variables states [1-5 . A well- 
known representative for the class of optical entangled 
states is the two-mode squeezed vacuum [6] which has 
been used as a resource for continuous variable quan- 
tum computation, cryptography and teleportation exper- 
iments |7fi9^. 

In the microwave domain, state tomography experi- 
ments have recently been realized for single itinerant field 
modes p!QHT2] . In addition, multi-mode state reconstruc- 
tion for intra-cavity fields of two spatially separated cav- 
ities has been demonstrated [13 . Here we present het- 
erodyne state tomography techniques that allow for the 
state reconstruction of two itinerant modes. We apply 
these techniques to reconstruct the Wigner function of 
a two-mode correlated state generated in a Josephson 
parametric amplifier. 

The general interest in parametric amplification [Mj 
D3] has steadily grown in the recent past jTBHTF because 
it allows to investigate the quantum properties of mi- 
crowave radiation [19], superconducting qubits [20] and 
nanomechanical oscillators [ST at high signal-to-noise- 
ratios. Noise squeezing at microwave frequencies has 
been demonstrated for these amplifiers in the degenerate 
as well as in the non-degenerate case of parametric 
amplification [22 . 

In this letter we present state tomography measure- 
ments for two output field modes of a parametric am- 
plifier. The two modes emitted from a single broadband 
resonator mode are separated in frequency space [23[ |24] . 
Complementary, one can analyze squeezing correlations 
between two spatially separated modes by sending two 
squeezed states trough a balanced beamsplitter [11 or 
by using two parametrically coupled resonators with dif- 
ferent output ports [22^. Here, we first describe the mea- 
surement setup, discuss the device parameters, and char- 
acterize the system as a phase insensitive amplifier. We 



then measure the covariance matrix and reconstruct the 
Gaussian Wigner function of the four quadrature compo- 
nents [7]. 

Our parametric amplifier [Fig. [l] is realized as a A/4 
transmission line resonator terminated by a supercon- 
ducting quantum interference device (SQUID) [^. The 
fundamental mode of the resonator is in good approxima- 
tion described as a nonlinear oscillator with an effective 
Hamiltonian 

K 

H = hhj^o) a ^ h—a^ aa^ (1) 

where the non-linearity is provided by the SQUID. The 
Kerr constant K as well as the resonator frequency uj^ /2t: 
are functions of the effective Josephson energy of the 
SQUID Thus, by changing the flux bias (j) through 
the SQUID loop the resonator frequency and the Kerr 
constant can be tuned. A measurement of the flux de- 
pendent resonance frequency allows to extract both the 
Josephson energy and the Kerr constant. Fitting the 




FIG. 1: Schematic of the experimental setup. The parametric 
amplifier is realized as a pumped A/4 resonator with a SQUID 
at the end. The input and and output modes are separated 
using a circulator. The output signal is amplified by a cold 
HEMT amplifier at the 4K stage and split into two channels 
at room temperature. Both channels have their individual 
local oscillator allowing for a simultaneous detection of two 
distinct modes in frequency space. 
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FIG. 2: (a) Measured (dots) and calculated (solid line) res- 
onance frequency cJr/Svr vs. magnetic flux 0. (b) Calculated 
real (blue line) and imaginary (red line) part of the reflection 
coefficient V in the linear regime for magnetic flux = and 
measured data (dots), (c) Measured maximum gain Co as a 
function of pump frequency ujp/2n and pump power Pp. A 
horizontal cut through the data where Pp ^ Pp,crit is shown in 
white, (d) Measured gain Ca (red) and power spectral den- 
sity 5'a (blue) for a fixed pump tone as a function of signal 
detuning A. The frequency filter functions /i(A) and /2(A), 
which define the modes hi and ^2, are shown on a logarithmic 
scale (arb. units) as the lines enclosing the shaded areas. 

data shown in Fig. [2|a) to a model described in Ref. [26] 
leads to the approximate values £^j,max /h ^ 6.1 THz and 
^/^r,max ^ —2.8 X 10~^ at the maximum resonance fre- 
quency where a;r,max/27r ^ 6.9 GHz. To further char- 
acterize the resonator we have measured the reflection 
coefficient F in the linear regime [Fig. [2|b)]. From the 
real and imaginary part of F we can extract the coupling 
of the resonator to the transmission line kI2ix 25 MHz 
which dominates over the internal loss yi/2iT ^ 2 MHz. 

If a nonlinear oscillator is pumped close to its bifurca- 
tion point with a coherent pump tone at frequency cjp/27r 
the relation between input and output field modes 6in(A) 
and 6out(A) can be written as [27 

6out(A) = AA6i„(A) + BA&U-A) . (2) 

Here A is the detuning from the pump frequency Co'p/27r, 
where A > corresponds to frequency components at the 
upper sideband of the pump whereas A < corresponds 
to the lower sideband. The frequency dependent coeffi- 
cients A A and 5a fulfill the relation |AaP — I^Ap = 1- 
Thus Eq. ([2| describes a minimal form quantum linear 
amplifier with gain Ga = I^Ap [15 . When a signal is 
applied at one sideband the frequency components of the 
other sideband are usually called idler modes. The de- 
vice acts as a phase-insensitive amplifier when only the 
signal modes are analyzed at its output [28 . 

We have characterized our device as such an ampli- 
fier. To measure the maximum gain Go as a function of 



u;p/27r and pump power Pp we apply an additional weak 
coherent tone with a small detuning of A/27r = 2.5 kHz 
to the input port. We measure the reflected amplitude 
at this frequency and compare the result with the one 
that we obtain when the pump tone is turned off. The 
absolute square of this ratio is the gain Go, see Fig.[2jc). 
We identify a critical point at cjp^crit/^Tr ^ 6.882 GHz 
and Pp,crit ^ — 80.6dBm, where Go takes its largest 
value [29]. For pump powers below Pp,crit we are in the 
stable amplifler regime. A decrease in pump power leads 
to smaller gain but to larger amplifler bandwidth B. The 
gain-bandwidth product remains constant according to 
the relation ^/GqB cx K: [27 , which we have verifled ex- 
perimentally. 

For the following measurements we have flxed the co- 
herent pump at ujp/27T = 6.8834 GHz and Pp ^ -80.8 
dBm. The signal frequency dependence of the gain Ga 
for this pump tone is shown in Fig. [2jd) together with 
a measurement of the power spectral density Sa- The 
power spectral density can be decomposed into two con- 
tributions 

5a<^(A-A') = {6L(A)6out(A'))+A^noise<5(A-A') . (3) 

The flrst term on the right hand side stands for the noise 
at the output of the parametric amplifler. The second 
term describes the white noise contribution of the HEMT 
amplifler. Since the power spectral density at the para- 
metric amplifler output can be interpreted as the am- 
plifled vacuum noise, the curves for Ga and Sa are al- 
most identical up to the constant offset A^noise = 69. In 
the large gain limit both curves are well described by a 
Lorentzian function [27 . 

Up to now we have characterized our device as a phase 
insensitive amplifler. If the input of the parametric am- 
plifler is in the vacuum state, which is approximately real- 
ized with a mean thermal photon number below n ~ 0.05 
as verifled in independent experiments [30 , the output 
at a single frequency component is the amplifled vacuum 
noise. However, if one analyzes the signal and idler modes 
individually the parametric amplifler transforms any in- 
put state into a state where the signal and idler modes 
show squeezing correlations. We observe these correla- 
tions for a vacuum input state by pumping the resonator 
with the same coherent tone as before and recording the 
quadrature amplitudes in both modes individually. 

The resonator output is further amplifled with a 
HEMT amplifler and split into two channels. To sep- 
arate signal and idler frequency components from each 
other [Fig. [l] we effectively down-convert the microwaves 
in both channels by mixing them with local oscillator 
tones at frequencies cjLOi/27r and cjL02/27r set 5 MHz 
above and below the pump frequency, respectively. The 
voltages in both channels are digitized every 10 ns with 
an analog to digital converter (ADC). Using a fleld pro- 
grammable gate array (FPGA) the data is digitally fll- 
tered with a symmetric fllter function /(A). As a result 
the two detection channels capture information about 
frequency components in the windows /i(A) = /(A — 
27r X 5MHz) and /2(A) = /(A + 27r x 5MHz), respec- 
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FIG. 3: Measured histograms for the quadrature pairs 
{Xi,X2} and {Pi,P2} with (a) pump tone on and (b) off. 
The gray elhpses indicate Unes of constant probabihty den- 
sity. 

tively [Fig. |2|d)]. The filter is designed such that both 
/i(0) = /2(0) = 0, filtering out the coherent pump tone. 

As a result, the complex numbers Si and S2 that we 
extract after the digital data processing correspond to 
measurement results of the complex valued operators [31] 

^l,2=^l,2+/^l,2=^l,2+^A,2 (4) 

where 

/oo 
dA/i,2(A)6out(A) . (5) 
-00 

hi and 62 describe a pair of signal and idler modes at 
the parametric amplifier output. The modes hi and /12 
are in thermal states with approximately the same mean 
photon number A^noise due to the noise added by the de- 
tection chain. 

Taking the filter function into account we get an op- 
erator relation similar to Eq. ([2| for 61, 62 and their 
counterparts at the input of the parametric amplifier. 
Such an operator identity is equivalent to the two-mode 
squeezing transformation described by the unitary oper- 
ator U{r) = ex.p[r{bib2 — ^1^2)] i^^.- relative phase 
between the two local oscillators has been chosen such 
that the squeezing parameter r is real and given by 

/oo 
dA|/i(A)|2GA . (6) 
-00 

From our gain measurement data and the designed fil- 
ter function we obtain the numerical value r ^ 1.75. 
A parametric amplifier device thus ideally transforms a 
vacuum input state into a two-mode squeezed vacuum 
state /7(r)|00) which shows strong correlations among 
the 4 quadrature components G ^2, ^2} of the 

two modes. Those are defined by 61^2 = ^1,2 + ipi,2- 



The relative "position" xi — X2 and the total "momen- 
tum" pi -\- p2 are squeezed according to {{xi — ^2)^) = 
{{pi + ^2)^) = e~^^/2, whereas each component itself 
is amplified {(^f) = cosh(2r)/4. Since the two- mode 
squeezed state belongs to the class of Gaussian states, 
its Wigner function can be written as a multivariate nor- 
mal distribution [7] 

with the vector of quadrature components a = 
^2,^2) and the quadrature covariance matrix V 
with elements Vij = {^iij + ijii)/2 [7]. The four- 
dimensional Wigner function of the two-mode phase 
space distribution is thus fully determined by the 4x4 
covariance matrix V. The elements of this matrix de- 
scribe the joint statistics of the amplitude fluctuations of 
the two modes. 

To determine these elements we detect the four quadra- 
ture components as explained above and store the re- 
sults in two-dimensional histograms for the six possible 
pairs {Xi,Pi}, {X2,P2}, {Xi,P2}, {X2,Pi}, {Xi,X2} 
and {Pi,P2}. For each pair we first acquire a refer- 
ence histogram with the pump turned off, which charac- 
terizes the quadrature distribution of the effective noise 
modes hi and /12 and a second histogram with the pump 
turned on [Fig. [3]. The single mode histograms {Xi, Pi}, 
{X2,P2} (data not shown) remain circular symmetric 
whereas the cross-histograms {Xi,X2} and {Pi,P2} are 
squeezed along the diagonal axis. Because of the rela- 
tively large added amplifier noise the differences between 
measurements with pump on or off are small but never- 
theless directly visible in the bare histograms. To scale 
X and P axes of the histograms we have used the infor- 
mation about the added amplifier noise extracted from 
the measurement of the power spectral density ^a- The 
axes are scaled such that the histogram variances in the 
"pump off' measurements are equal to the noise offset 
Xnoise- Furthermore, we have corrected for the small dif- 
ference of gain in the two detection channels. 

The measured data is further analyzed by separating 
the contributions of modes hi and /12 from those of modes 
bi and 62- We calculate all possible expectation values 
(X^Xj), (XiPj) from the 12 measured histograms. Tak- 
ing into account that modes hi and bj are uncorrelated 
and using Eq. 0, we can determine all second order ex- 
pectation values {^i^j) fiO\ fT2], which are summarized 
in the covariance matrix shown in Fig. |4ja). The diago- 
nal elements in the matrix express the amplified individ- 
ual quadrature fluctuations in both modes, whereas the 
non- vanishing off-diagonal elements describe the squeez- 
ing correlations between the two modes. If we fit our 
measurement results to the ideal values for which the di- 
agonal terms are cosh(2r)/4 and the four non- vanishing 
off diagonal terms ±sinh(2r)/4, we find r ~ 1.78, which 
is close to the value that we have calculated from the in- 
dependent gain measurement data and the designed filter 
function (r ~ 1.75). 
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FIG. 4: (a) Measured quadrature covariance matrix V of 
the two-mode correlated state which allows to reconstruct the 
four-dimensional Wigner function, (b) and (c) show projec- 
tions of this function onto two-dimensional subspaces. (d) For 
comparison, we show the same projection of the theoretical 
Wigner function of an ideal two-mode squeezed vacuum state. 

We note that the measured diagonal elements are ap- 
proximately 3% larger and the off-diagonal elements the 
same amount smaller than the ideally expected ones. The 
output of the parametric amplifier in modes bi and 62 is 
thus not perfectly described by the probability distribu- 
tion of a pure two-mode squeezed vacuum state but by 
its convolution with a thermal distribution. This thermal 
noise adds to the individual quadrature fluctuations and 
reduces the effective squeezing. Finite internal loss is one 
of the potential sources of uncorrelated noise [27 . 

We evaluate Eq. ^ to reconstruct the four- 
dimensional Wigner function W{a) for the two modes. In 
Fig. [3|b) and (c) we show projections of W{a) on two- 
dimensional subspaces. The {xi,pi}-projection, which 
describes the individual state of mode 61, is thermal. 
Compared to the vacuum Wigner function (not shown) 
it has a larger variance, indicating the amplified vacuum 
noise. The {xi, X2}-projection is squeezed along the di- 
agonal axis, visualizing the correlations between the two 
modes. Due to the additional uncorrelated noise the re- 
constructed distribution is broader than it would be in 
the case of an ideal squeezed two-mode vacuum state (see 
Fig. [sj^d)) as explained above. 

In summary, we have reconstructed the full quadrature 
covariance matrix of a two-mode correlated state to de- 
termine its four-dimensional Gaussian Wigner function. 
In this context, we have fully characterized a Josephson 
parametric amplifier device. Our experiments show that 
linear photon detection allows for efficient state tomog- 
raphy measurements of two radiation field modes. We 
have observed a squeezing parameter of r ^ 1.78 and be- 
lieve that this value can easily be increased by building a 
parametric amplifier with larger bandwidth. This should 
improve the ratio between amplifier and detection band- 
width which furthermore allows to use the device as a low 
noise amplifier in other circuit QED experiments [20]. 



In addition, combining parametric amplifier devices in 
networks with beamsplitters and superconducting qubits 
could allow for future continuous variable quantum com- 
putation with propagating microwave photons 0. 
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